Magnetic interaction between GaMnAs layers via spin-polarized 
quasi-two-dimensional hole gas: Monte Carlo simulation 
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The magnetic order resulting from an indirect exchange between magnetic moments in the metal- 
lic phase of a Gai-2,Mn3;As bilayered structure is studied via Monte Carlo simulation. The coupling 
mechanism involves a perturbative calculation in second order of the interaction between the mag- 
netic moments and carriers (holes). We take into account a possible polarization of the hole gas 
establishing, thus, self-consistency between the magnetic order and the electronic structure. It leads 
to a ferromagnetic order even in the case of thin layers. This fact is analyzed in terms of the inter- 
and intra-layer interactions. 
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The possibility of having a diluted magnetic semi- 
conductor (DMS) based on GaAs opens a wide range 
of potential applications as, for instance, in integrated 
magneto-optoelectronic devices In Gai_j;Mna;As al- 
loys substitutional Mn acts as an acceptor (it binds one 
hole), and at the same time it carries a localized mag- 
netic moment, due to its five electrons in the 3d shell. 
For X — 0.053, the alloy is a metallic ferromagnet the 
Curie- Weiss temperature is 11 OK, and the free hole con- 
centration is near 10^*^0771"^. The ferromagnetic order in 
the metallic phase is understood, at present, as resulting 
from the indirect exchange between the Mn^+ ions due 
to the local spin polarization in the hole gas. This ex- 
planation implies the spin coherence length to be larger 
than the average distance of the localized moments. 

In this work we extended a confinement- adapted 
RKKY mechanism to study the magnetic order re- 
sulting of the indirect exchange between magnetic mo- 
ments in a GaAs/Gai_a;Mna;As nanostructure with two 
DMS layers. A Monte Carlo simulation is performed to 
determine the resulting magnetic phases. The interac- 
tion potential between a Fermi gas and a set of localized 
magnetic moments is well described by the Hund-type 
exchange potential: 
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where Si is the localized spin of the Mn ion at position 
which will be treated as a classical variable, and s{r) 
is the spin operator of the fermion at r; / is the sp — d 
interaction In terms of the hole field operators, the 
Hund Hamiltonian is 



Instead of free fermions in a 3-D space, the electrons and 
holes in a semiconductor hetcrostructure are confined in 
the growth direction, assumed to be the z-axis, due to 
the mismatch of the conduction and valence band edges. 
The total Hamiltonian, H — Hq+Hc^ includes the kinetic 
part, the confinement potential and the Hartree as well 
as exchange and correlation terms in Hq |^ . 

Neglecting the scattering by impurities, holes are free 
particles in the plane perpendicular to that growth di- 
rection, i.e., in the plane parallel to the semiconductor 
interfaces. Their field operator can be written as 
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where A is the normalization area, fc is a wave vector in 
the plane {x,y), rj^ is the spin tensor for the polarization 
CT, <f>n,cr{z) is the envelope function which describes the 
motion of the fermion in the z-direction, and j: ^ is the 

fermion annihilation operator for the state {n,k,cr). p 
represents a vector in the 2-D coordinates plane {x,y). 

The confined RKKY indirect exchange is a second or- 
der perturbative treatment This approach can be 
extended to start the perturbative calculation from spin- 
polarized states, going, in one sense, beyond the sec- 
ond order of perturbation. In that case, we do not ob- 
tain a simple scalar product of localized dipole moments 
in a Heisenberg-like Hamiltonian, since the polarization 
breaks the rotational symmetry establishing a preferen- 
tial direction, that of the average magnetization. After 
a lengthy calculation, however, it can be shown that the 
interaction term can be written as an effective Hamilto- 
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FIG. 1: Spin-polarized charge density distributions for two 
Gao.95Mno.05As layers of lOA , separated by a 20A width 
GaAs layer, and hole concentration p = 1 x lO^^cm"^. Solid 
line for total charge density, dotted line for spin t (parallel), 
dashed line for spin I (anti-parallel). The polarization density 
in dotted-dashed line. 

{cjl + cj:])is:s^ + sfsy) (4) 

where x, y and z (the later assumed to be the direction 
of magnetization) are spin coordinates, and i and j are 
indices of the magnetic moment. The sp — d interaction, 
is generally written as I/vq = NqP, where vq is the effec- 
tive volume of the Mn ion {uq/Ax in the present case), 
A'o is the density of atoms of Ga, and /3 is the sp — d ex- 
change. For GaMnAs, NqP = -1.2eV §]. The exchange 
coefficients C,^'' are expressed in terms of the real space 
Fourier transform of the Lindhard function: 

= - E E E(^)''^"'(^»)<^»(^^)x 

x</.:(z,)0„,(z,)x"'"'(R.,) (5) 

with 

X"'"'(R,;,) = X^exp[-*g. R.,]x"'"'(g) (6) 

k 

The Lindhard function is defined by 0: 

A.""'" (g) = " n',k+q' 

We have focused our attention on bilayered systems, 
consisting of two GaMnAs layers inside an otherwise non- 
magnetic GaAs matrix. Fig.(l) shows the density of 
holes, charge density for each spin orientation, and the 
polarization density for the case where two lOA GaMnAs 
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FIG. 2: Magnetization as a function of temperature for two 
Gao.95Mno.05As layers of lOA width and hole concentration 
p = 1 X lO^^cm"^, separated by a GaAs layer whose width in 
A is indicate by arrows. 

layers separated by 20A have their magnetic moments 
fully aligned. It results from a self-consistent calcula- 
tion Q of the eigenstates of with x — 5% and the 
hole concentration is 10 x 10^°cm~^. The interaction, 
given by Eq.(^, is assumed to be effective within a cut- 
off radius, which we have taken as i?c = 2.5a, where a 
is the fee lattice parameter of GaAs. The calculation 
is performed in a finite box, whose axes are parallel to 
the [100] directions. Its dimensions are = Ly, and 
Lz — Na/2, and N is the number of DMS monolayers 
(ML). Periodic boundary conditions are imposed in the 
(x, y) plane. The lateral dimensions are adjusted in such 
a way that the total number of spins ranges from 
3500 to 4100, for all samples with different L^. The ini- 
tial spin orientations are randomly assigned. At a given 
temperature, the energy of the system due to the in- 
teraction is calculated, and the equilibrium state for a 
given temperature is sought by changing the individual 
spin orientation according to the Metropolis algorithm. 
A slow cooling stepwise process is accomplished making 
sure that the thermal equilibrium is reached at every tem- 
perature. The resulting spin configuration is taken as the 
starting configuration for the next step at a lower tem- 
perature. For every temperature, the average magnetiza- 
tion < M > is calculated. Whenever the magnetization 
reaches thresholds of multiples of tenths of 5?i/2, the elec- 
tronic structure is recalculated assuming a homogeneous 
magnetization in the DMS layers, and the G^^ coeffi- 
cients are also recalculated according to Eq. (a). In a 
sense, this approach establishes a self-consistency in the 
calculation of the magnetic order and in the electronic 
structure. 

We have tested the influence of the separation of the 
two DMS layers. It is known that a ferromagnetic order 
can not occur in a single DMS layer below a certain layer 
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FIG. 3: Same as above for p = 2 x 10^° cm"^ 

width |3|. However, the magnetic order in a DMS bilayer 
responds to the interaction occurring both intra- and 
inter- layer. Fig. (2) shows the magnetization as a func- 
tion of temperature, obtained for DMS thickness of 10 A 
in each layer, with a hole concentration of p = 1 x 10^° 
cm^'^, for different layer separation. We have chosen an 
interaction cutoff radius of 5 monolayers, that is greater 



them the mean Mn distance, and large enough to take ac- 
count for anti-ferromagnetic interactions. Increasing the 
inter-layer separation, the Curie temperature decreases 
as a direct consequence of the fact that the inter-layer 
interaction decreases with the separation. Fig. (3) shows 
results for a similar calculation, but p = 2 x 10^° cm^'^. 
We observe a rising on the Curie temperature. Also, in- 
creasing the carrier concentration implies in increasing 
the Fermi wave vector what introduces the possibility 
of anti-ferromagnetic interactions within the chosen cut- 
off radius. This produces a partial magnetization in the 
curves, associated with the occurrence of canted phases 

Our results show that a ferromagnetic order with sig- 
nificant Curie temperature occurs in bilayered structure, 
even in the case of thin DMS layers. For the sake of com- 
pleteness we calculated also some samples with two DMS 
layers of 20 A each. On these systems we observed the 
same behavior described above, but with higher Curie 
temperatures as expected for thicker layers. 
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